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Abstract 

The Dirac equation for an electron in two spatial dimensions in the Coulomb and 
homogeneous magnetic fields is a physical example of quasi-exactly solvable systems. 
This model, however, does not belong to the classes based on the algebra sl{2) which 
underlies most one-dimensional and effectively one-dimensional quasi-exactly solvable 
systems. In this paper we demonstrate that the quasi-exactly solvable differential 
equation possesses a hidden osp{2, 2) superalgebra. 
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1 Introduction. 



It is well known that exactly solvable models play an important role in many fields of 
physics. However, exactly-solvable systems are very rare in physics. Recently, in quantum 
mechanics a new type of spectral problems, which is intermediate to exactly solvable ones 
and non solvable ones, has been found [[l]-[7]]. For this new class of spectral problems, 
the so-called quasi-exactly solvable (QES) models, it is possible to determine algebraically 
a part of the spectrum but not the whole spectrum. Such quasi-exact solvability is usually 
due to some hidden Lie-algebraic structures. More precisely, a QES Hamiltonian can be 
reduced to a quadratic combination of the generators of a Lie group with finite-dimensional 
representations. 

The first physical example of QES model is the system of two electrons moving in an 
external oscillator potential discussed in [8, 9]. Later, several physical QES models were 
discovered, which include the two-dimensional Schrodinger [10], the Klein-Gordon [11], and 
the Dirac equation [12, 13] of an electron moving in an attractive/repulsive Coulomb field and 
a homogeneous magnetic field. More recently, the Pauli and the Dirac equation minimally 
coupled to magnetic fields [14], and Dirac equation of neutral particles with non-minimal 
electromagnetic couplings [15] were also shown to be QES. 

It turns out that the system of two electrons moving in an external oscillator potential 
and those of an electron moving in an Coulomb field and a homogeneous magnetic field 
mentioned above, with the exception of the Dirac case, share the same underlying structure 
that made them QES. In [16] it was demonstrated that these systems are governed essentially 
by the same basic equation, which is QES owing to the existence of a hidden s/(2) algebraic 
structure. This algebraic structure was first reahzed by Turbiner for the case of two electrons 
in an oscillator potential [17]. For the Dirac case, on the other hand, it had been shown the 
quasi-exact solvability of this system is not related to the sl{2) algebra [13]. At that time, 
however, it was not known to us what algebraic structure the Dirac system possesses. In 
this paper we would like to show that the relevant symmetry of this system is the osp{2, 2) 
superalgebra. 
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2 The Dirac Equation 

In 2 + 1 dimension, the Dirac equation for an electron minimally coupled to an external 

electromagnetic field has the form (we set c = h = 1) 

{idt-HD)^{t,r)^0, (1) 

where 

Hd = <JiP2 - (J2P1 + (73m - eA° (2) 

is the Dirac Hamiltonian, ak {k — 1,2,3) are the Pauli matrices, Pk — —idk + eA^ is the 
operator of generalized momentum of the electron, the vector potential of the external 
electromagnetic field, m the rest mass of the electron, and — e (e > 0) is its electric charge. 
The Dirac wave function ^(i, r) is a two-component function. For an external Coulomb field 
and a constant homogeneous magnetic field B > along the z direction, the potential 
are given by 

A°(r)=Ze/r (e>0). A, ^ -By/2, Ay^Bx/2 (3) 
in the symmetric gauge. The wave function is assumed to have the form 

x) = ^ eM-^Et)^Pi{r, if) , (4) 
where E is the energy of the electron, and 



(5) 



with integral number I. The function ipi{r,(fi) is an eigenfunction of the conserved total 
angular momentum Jz = Lz + Sz = —id /dip + (73/2 with eigenvalue j = / + 1/2. It should 
be reminded that / is not a good quantum number. Only the eigenvalues j of the conserved 
total angular momentum are physically meaningful. 

Putting Eq.(4) and (5) into (1), and taking into account of the equations 



we obtain 



dF 

dr 
dG 



'l + l eBr\ ^ / Za\^ ^ 



+ 



+ \ eBr\ ^ 



E-m + F = , 

r / 



(7) 
(8) 



dr \ r 2 ^ 

where a = e'^ — 1/137 is the fine structure constant. In a strong magnetic field the asymp- 
totic solutions of F{r) and G{r) have the forms exp(— at large r, and r'^ with 
7 = ^{l + l/2y — {Zay for small r. One must have Za < 1/2, otherwise the wave function 
will oscillate as r — > when I — and I — —1. With these asymptotic factors, we now write 



F{r) = r^exp(-eSrV4) Q{r), G{r) = exp(-eSrV4) P{r) . 



(9) 



Substituting Eq.(9) into Eq.(7) and (8) and eliminating P(r) from the coupled equations, 
we obtain 



dr' 



+ 



27 



Za/r2 

QBr H 

r E + m + Za/r 



dr 



+ 



2EZa I + h 



+ ^-5-e^(r + i) 



Za/r^ 



eBr 



I + 1/2' 



Q{r) = , 



(10) 



E + m + Za/r r 

where F = / + 1/2 + 7. Once Q{r) is solved, the form of P{r) is obtainable from Eqs.(7) and 
(9). If we let X = r/Zfi, = 1/VeB, Eq.(lO) becomes 



d' 



+ 



dx^ 
HE' 



27 Za 

— ■ — X -\ 

X x{{E + m)lBX + Za) 



d_ 

dx 



Za{l + 1/2 - 7) 



X 



x^ 



Za 



x'^[{E + m)lBX + Za\ (E + m)lBX + Za 
Eq.(ll) can be rewritten as 



Qix) = 



(11) 



dx'^ 



+ 



2/? 



X 



1 



X 



X + Xo 



d be 

— + e + 

dx X X + Xq 



Q{x) = . 



(12) 



Here /3 = 7 + 1/2, Xq = Za/[{E + m)/^], e = {E^ - m^l - {T + 1), b = bo + L/xo, 
bo — 2EZalB, L = {I + 1/2 — j), and c = Xq + L/xq. The energy E is determined once 



the values of e and xq are known. The corresponding value of the magnetic field B is then 
obtainable from the expression Is = Za/[{E + m)xo]. Solution of e and xq is achieved in 
[13] by means of the Bethe ansatz equations. There it was shown that Eq.(12) is QES when 

e is a non-negative integer, i.e. e — n, n — 0,1,2 In this case, the function Q{x) is a 

polynomial of degree n. 

If we eliminate Q{x) from Eqs.(7) and (8) instead, we will obtain a second differential 
equation of P{x): 



213 1 

X 7 

X X -\- X r: 



, + e' + \ P{x) = , (13) 

dx X X -\- Xq\ 



with = Za/[{E — m)/^], e' — e -\- 1, h' — ho -\- d , and c' — —V/x'q. Other parameters are 
as defined previously. It is obvious that Eq.(13) is in the same form as Eq.(12), and hence 
is also QES. In [13] it was shown that Eq.(12) and (13) gave the same QES spectrum when 

e — n , e' = n+l, n = 0,1,2,... , 
b' -c' = b-c + xo , (14) 

The result that e = n and e' = n+1 implies that the degree of the polynomial P{x) is of one 
order higher than that of Q{x), i.e. when e — n, Q{x) and P{x) are of degree n and n + 1, 
respectively. 

In [13] it was also shown that the quasi-exact solvability of Eqs.(12) and (13) is not 
due to the sl{2) Lie-algebra, which is responsible for the quasi-exact solvability of most one- 
dimensional and effectively one-dimensional systems. We now want to show that the relevant 
algebra is indeed osp{2, 2). In what follows, we shall present a differential representation of 
this algebra, and then use it to demonstrate how osp{2, 2) underlies the QES structure of 
our Dirac system. 
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3 A Differential Representation of osp(2, 2) 



The superalgebra osp(2, 2) is characterized by four bosonic generators T^'^, J and four 
fermionic generators (^1,2, Qi,2- These generators satisfy the commutation and anti-commutation 
relations [4, 7] 

[r°, T^] = ± , [T+, T-] - -2r° , [J, T"] = , a = 0, +, - , 
{Qi,Q2}^-T- , {Q2,Qi} = T+ , 



^({ai,Qi} + {a2,g2}) = J , l{{Qi,Qi}-{Q2,Q2})^T' , 



[Qi,T+]=Q2 , [Q2,T^]=0 , [Qi,T-] = , [Q,,T-] = -Q, , (15) 



[gi,2,T°] = ±^Ql,2 , [Qi,2,T'] = tIQi,2 , 



1 - 1 - 

[Ql,2, J] — —-^Ql,2 , [Ql,2, J] = 2^1,2 ■ 



A differential representation of the algebra osp{2, 2) can be reahzed by the following 2x2 
differential-matrix operators: 







x^dx — {n+ l)x 



rpO 



rd —- 



T: 



'^4 ^ 

y 4 y 



^ _n±2 n ^ 

2 " 

n —Mil 

^ 2 / 



(16) 
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Qi = 



1 



^2 = 



X 








(n + 1) 







02 = 



— Or 







Here x is a real variable, = d/dx, and n is a real number. It is easily checked that 
these matrices satisfy the osp(2, 2) (anti-) commutator relations Eq.(15). For non-negative 
integer n , there exists for the osp{2, 2) algebra a (2n + 3)-dimensional representation, with 
an invariant subspace P^^i consisting of two-component functions of the form 



ip{x) 



qn[x) 

Pn+l(x) j 

where and are polynomials of degree n and n + 1, respectively. 



V 



(17) 



4 Hidden Lie- Algebraic Structure of The Dirac Equa- 
tion 

We now show that the osp(2, 2) super algebra is indeed the hidden algebraic structure un- 
derlying the quasi-exact solvability of Eq.(12) when e is a non-negative integer n > 0. First 
we rewrite Eq.(12) as 

Tq{x) Qn{x) = , 

Tq{x) = {x'^ + Xqx)— + (-x^ - xox^ + 2/3 X + 2/3 xo) — 

-|- nx'^ -\- {nxQ -\-h — c)x -\- bxo . (18) 

We have put e — n, and write Q{x) as Qn{x) to indicate that it is a polynomial of degree n. 
This equation may be cast into the matrix form 



Qn{x) 


\ 


J 


"1 


Qn{x) 


] 




(19) 










= . 




/ 




% 




) 
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The operator Tq{x) is the 2x2 matrix operator in Eq.(19). Written in this form, the 
super algebraic structure hidden in it can be clearly exhibited. The 2x2 matrix operator 
Tq turns out to be expressible as a linear combination of the osp{2, 2) generators in Eq.(16). 
The two-component wave function in Eq.(19) is simply an element in the invariant subspace 
Pn+i °f algebra with the lower component Pn+i{x) = (the form of Pn+i{x) is immaterial 
here, as Tq annihilates the lower component of any element in the subspace P^+i) ■ We 
shall demonstrate this below. 

First let us express the various terms in Tq in terms of the generators in Eq.(16). After 
some algebras, we obtain the following results: 



( o^ 



x'^dl 



Q^T-T- 



( o^ 



xdl 



x^dx — nx'^ 



( 0^ 



x^dx 



2{Q{rl Q,T+) 



^ o^ 



nx 



f 0^ 



xd^ 



QiT- 





d^ 



With these expressions, the operator Tq can then be written as 



Tq = {2Q2T^T- - QiT^T-) + X0Q2T-T- - 

-xoi2Q2T^ - QiT;t) + 2f3iQ2T- + xoQiT') 
+2xo{Q2T^-QiT+) + {b-c)Q2 + bxoQi . 



(20) 



We have therefore succeeded in expressing Tq as a hnear combination of the generators of 
osp{2, 2) with finite dimensional subspace. The underlying algebraic structure responsible 
for the quasi-exact solvability of Eq.(12) is thus demonstrated. 

The underlying algebraic structure responsible for the quasi-exact solvability of Eq.(13) 
can be obtained in the same way. In fact, as mentioned in Sect. 2, Eq.(13) is in the same 
form as Eq.(12), only with the parameters x, b, c and n being replaced by x', b', c' and 
n+1, respectively. Hence, we immediately see that the corresponding operator Tp will be of 
exactly the same form as Tg with the same replacements of the corresponding parameters. 

5 Summary 

In this paper we have unveiled the underlying symmetry responsible for the quasi-exact 
solvability of planar Dirac electron in a Coulomb and a magnetic field. The relevant second 
order differential operator acting on any one component of the two-component wave function 
was recast in a 2 x 2 matrix form. This 2x2 differential-matrix operator was then shown to 
be expressible as a linear combination of the generators of the super algebra osp(2,2), thus 
exhibiting the algebraic structure of the QES Dirac system. With this result, all the algebraic 
structures making the systems of planar charged particles in Coulomb and magnetic fields 
QES have been identified. 
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